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A derivative nonlinear Schrodinger model is shown to support localized A-body bound states 
for several ranges (called bands) of the coupling constant r\. The ranges of r\ within each band can 
be completely determined using number theoretic concepts such as Farey sequences and continued 
fractions. For TV > 3, the iV-body bound states can have both positive and negative momentum. 
For r\ > 0, bound states with positive momentum have positive binding energy, while states with 
negative momentum have negative binding energy. 
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Integrable quantum models in 1+1 dimensions which support bound states have been studied extensively for many 
years [1-8]. For an integrable Hamiltonian, the coordinate Bethe ansatz can yield the exact eigenfunctions. If an 
eigenfunction decays exponentially fast when any of the interparticle distances tends towards infinity (keeping the 
center of mass coordinate fixed), we call such a localized square-integrable eigenfunction a bound state. Bound states 
^. , of quantum integrable models are usually found to have positive binding energy [1-5] . 

■^j- ■ In this paper, we will study the quantum bound states of an integrable derivative nonlinear Schrodinger (DNLS) 

model [6,7]. Classical and quantum versions of the DNLS model have found applications in many areas of physics 

~H , like circularly polarized nonlinear Alfven waves in a plasma [9], quantum properties of solitons in optical fibers [10], 

and some chiral Luttinger liquids [11]. The classical DNLS model is known to have solitons with momenta in only 

one direction [12,13]. 

Using the coordinate Bethe ansatz, it had been found earlier that quantum A^-body bound states exist for the DNLS 
model provided that the interaction parameter r\ (defined in Eq. (1) below) lies in the range < |7y| < tan(7r/A r ). It 
was also observed that, similar to the classical case, such A^-body bound states can have only positive values of P/rj, 
where P is the momentum [6]. However, it was found recently that bound states can exist with Pjr\ < provided 
j-~< that ta,n(n/N) < \r/\ < tan[7r/(A~ — 1)], and that these states have negative binding energy [8]. This naturally leads 
one to ask: are there other ranges of values of r\ for which quantum bound states exist, and, if they exist, what are 
their momenta and binding energies? 

In this paper, we will solve the problem of determining the complete ranges of values of 77 for which quantum 
A~-body bound states exist in the DNLS model, for all values of N. After presenting the conditions which are required 
for a quantum A^-body bound state to exist, we will use the idea of Farey sequences in number theory to show that 
there are certain ranges of rj, called bands, in which bound states exist. We find that the bound states appearing 
within each band can have both positive and negative values of Pji); these have positive and negative binding energies 
respectively. We will then use another concept from number theory, that of continued fractions, to address the inverse 
problem of finding the values of A~ for which A^-body bound states exist for a given value of 77. 

For A" particles, the Hamiltonian of the DNLS model is given by 

N 02 q f) 

"- = - "' E 51 + ^'" E «<* - «.) (hj + sj-) . CU 

j = l J l<m 

where we have set the particle mass ra = 1/2. Hn commutes with the momentum operator Pn = — ihy~) ,_ 1 d/dxj. 
We note that Hn remains invariant while Pn changes sign if we change the sign of 77 and transform all the x% —^ —Xi 



* E-mail: biru@tfieory.saha.ernet .in 
^ E-mail: tanaya@theory.saha.emet. in 
"'"E-mail: diptiman@cts.iisc.ernet.in 



at the same time; we call this the parity transformation. Hence it is sufficient to study the model for one particular 
sign of 77, say, 77 > 0. The cigenfunctions for 77 < can then be obtained by changing Xj — > — x,. 

Next, the coordinate space R N = {xi,X2, • • -Xn} is divided into various iV-dimensional sectors defined through 
inequalities like x^m < x w ( 2 ) < • • • < x w (j\n, where w(l),w(2), • • • ,w(iV) represents a permutation of the integers 
1, 2, • • • , N. Given the wave function in the fundamental sector, defined as x\ < x^ < ■ ■ ■ < X/v, the wave functions 
in all the other sectors can be found by Bose symmetry. In the fundamental sector, the Bcthe ansatz wave function 
takes the form 

ip = ^2 C w exp {i(fc w( i)Xi H hfc w(JV) x A r)}, (2) 

where fc„'s are all distinct wave numbers, the sum is over all permutations u) of the integers 1, 2, • • • , N, and C u are 
appropriate coefficients. The momentum and energy of this eigenfunction are given by 

N N 

p = nj2 k 3 1 and E = h2 J2tf ■ ( 3 ) 

In Ref. [6], it was shown that a localized bound state has a wave function consisting of only one plane wave in 
each sector. Namely, in the fundamental sector, the coefficients C u in (2) vanish for all w's except for the identity 
permutation. Further, the momenta fc„'s satisfy the following conditions: 

k n - k n +i + i V (k n + fc„+i) = , (4) 

for n= l,2,---,iV-l, 

N 

5> = 0, (5) 

3=1 

where qj denotes the imaginary part of kj , and 

N-l 

qi < , qi + q 2 < , • • • , ^2 q 3 - < . (6) 

Eqs. (5-6) imply that the wave function ip is square-integrable if one holds the center of mass coordinate X = J^ . Xi/N 
fixed, and integrates over the relative coordinates y r = x r +i — x r , where r = 1, 2, • • • , N — 1. In the fundamental 
sector, the integrals over the y r 's all run from to 00, and they are independent of each other. Using Eq. (5-6), one 
can express the probability density as 

JV-l r 

H 2 - exp[ J^ (E<U)Vr], (7) 

r=l j=l 

which is independent of X. Due to the conditions in (6), integration of this probability density over the y r 's gives a 
finite result. 

The conditions (4) and (5) imply that the k n 's must be complex numbers of the form 

k n = X e-^ 1 -^ , (8) 

where x is a real parameter, and <f> = tan -1 r\. Since < \i)\ < cxd, we assume that < \cf>\ < tt/2. Using Eqs. (3) and 
(8), the momentum and energy of this state are found to be 

^ sm(N<j)) , h 2 x 2 sm(2N<j)) 

p = hx —Hr ' and E = ■ (L\ ■ (9) 

If we define the mass of this state by the relation E = P 2 /(2M), we find that 

M = P^- . (10) 

2tan((/>) v ' 

We now impose the conditions (6) on the fc n 's. We find that all the following inequalities must be satisfied, 



X^f M(N~l)ct>] > (11) 

sm0 

for I = 1, 2, • • • , N - 1. For iV = 2, (11) is satisfied when lies in the range < < tt/2 (-tt/2 < < 0) if \ > ° 
(x < 0). Thus any nonzero value of <fr allows a 2-body bound state. From (9), we see that the ratio P/<j> > 0. 

We now consider the more interesting case with N > 3. Due to the parity symmetry of (1), we will henceforth 
assume that 4> > (i.e. n > 0). Eq. (11) can then be rewritten as 

X sin(^) sin[(N - 1)0] > (12) 

for I = 1, 2, • • • , N - 1. For x > 0, Eq. (12) implies 

cos[(N - 2l)(f>} > cos(N<f>) (13) 

for I = 1,2,- ■ ■ ,N — 1. Let us now consider a value of </> of the form 

0JV,n = ^ . ( 14 ) 

where n is an integer satisfying 1 < n < N/2. If n is odd, cos(N(f>N. n ) = — 1. We then find that all the inequalities 
in (13) are satisfied provided that N and n are relatively prime, i.e., if the greatest common divisor of TV and n is 1. 
Similarly, for x < 0, Eq. (12) takes the form 

cos[(N - 2l)<f>] < cos(N(t)) (15) 

for I = 1, 2, • • • , TV — 1. We find that all these inequalities are satisfied if n is even, and N and n are relatively prime. 

In short, all the inequalities in (12) are satisfied for (f> — 4>N.n, if and only if N and n are relatively prime (with n 
odd for x > 0, and n even for x < 0). By continuity, it follows that all the inequalities will hold in a neighborhood of 
4>N,n extending from a value <f>N,n,- to a value <j>N, n ,+- The region <f>N,n,- < 4> < <}>N,n,+ will be called the band B/v,n- 

For a given value of N, the number of bands in which bound states exist is equal to the number of integers n which 
are relatively prime to N and satisfy 1 < n < N/2. This is equal to half the number of integers which are relatively 
prime to N and satisfy 1 < n < N . The latter number is called Euler's 0-function $(iV) [14]. The number of bands 
is therefore equal to $(iV)/2 for -q > 0. 

We now have to determine the end points (j>N,n,— an d 4>N,n,+ of the band Sjv, n - One or more of the inequalities in 
(12) will be violated at the end points <f>N,n.± if 



PN,n,± 

where j± and l± are integers satisfying 



^ , (16) 



1 < l ± < N , and j ± < y (17) 

(since <p < 7r/2). Thus the end points of the band B/v,n are given by two rational numbers of the form j±/l± which 
lie closest to (and on either side of) the point 4>n n /ir — n/N. These can be found using the idea of Farey sequences 

[14]- 

For a positive integer N, the Farey sequence Fn is defined to be the set of all the fractions a/b in increasing order 
such that (i) < a < b < N, and (ii) a and b are relatively prime. For N > 2, if n/N is a fraction appearing 
somewhere in the sequence Fn, then it is known that the fractions a\/b\ and 02/62 appearing immediately to the left 
and to the right respectively of n/N satisfy 

o-i , a,2 < n , and a\ + 0,2 = n , 

h , b 2 < N , and b x + b 2 = N , 
nb\ — Na\ = 1 , and n&2 — Na,2 = — 1 , (18) 

and n, 61, 62 are relatively prime to AT. 

Using Eqs. (16) and (17), we now see that the end points of the band B^^ n are given by 

<PN,n,- = -, — , and <pN,n,+ = -; — , (19) 

01 b 2 
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FIG. 1. The values of <t>/n for which JV-body bound states exist for various values of N. 



where a\jb\ and CJ2/&2 are the fractions lying to the left and right of n/N in the Farey sequence F^. 

For N > 3, the lowest band is given by n = 1; by using Eq. (18), the range of this band is obtained as < <j>/ir < 
l/(iV — 1). For higher values of n, the end points of the band Bm,u (i.e., the integers a>i and bi) can be determined 
numerically by using the properties given in Eq. (18). Fig. 1 shows the ranges of values of <j> for which bound states 
exist for N = 2 to 20. 

Eq. (18) implies that the width of the right side of the band -B/v i?l from <pN,n to 4>N.n.+ is Tr/(Nb2), while the width 
of the left side from 0jv jrl] _ to 4>N,n is n/(Nb\). For later use, we note that each of these widths is larger than w/N 2 , 
since b\, b 2 < N. 

We now calculate the momentum and binding energy for the iV-body bound states in a particular band Bjy.n using 
Eq. (9). The form of the end points given in Eq. (19) shows that sin(iV</>) = at only one point in the band -B/v,n, 
namely, at <j> — <j>N,n- In the right part of the band (i.e., from <pN,n to (j>N, n ,+), the sign of sm(7V<^) is (—1)™. In the left 
part of the band (i.e., from 0jv in ,- to (j>N.n), the sign of sm(N<f>) is (— 1)" +1 . Since \ nas the same sign as (— 1)" +1 , 
the momentum given in Eq. (9) is positive in the left part of the band, negative in the right part of the band, and 
zero at </> = <pN,n- 

To calculate the binding energy, we consider a reference state in which the momentum P of the ./V-body bound 
state is equally distributed among TV single-particle scattering states. From Eqs. (3) and (9), the wave number 
associated with each of these single-particle states is found to be fco = xsin(./V0)/(7Vsin(/>). The total energy for the 
N single-particle scattering state is therefore given by 



E» 



h^Nki 



h\ 2 sm.\N<t>) 

N sin 2 6 



Subtracting E in (9) from E s in (20), we obtain the binding energy of the TV-body bound state as 

h 2 x 2 sin(iV</>) r sin(7V</>) cos(iV^) •> 



E B {<t>,N) 



sm< 



"I iVsir 



cos< 






(20) 



(21) 



Substituting N = 2 in Eq. (21), we obtain Eb((J>, 2) = 2h 2 x 2 sin 2 cj>. Thus Eb((J>, 2) > for any nonzero value of . 
Let us now consider the case N > 3. We can rewrite Eq. (21) in the form 



w 




(j)/7T 



FIG. 2. The binding energy Eb of the TV-body bound state as a function of (j)/ir for three different values of TV. 



E B (cj),N) 



h xsin(TV</)) 



TV sin tficos 

/(0,/V) = X [sin(/V0)cos0 



Ncos(N<j))sin(t)}. 



(22) 



AT-1 



On adding up all the inequalities given in (13) or (15), and using the identity X^=i cos[(TV— 2l)<fr] = sin[(TV— 1)<T)]/ sin cj>, 
we find that /((/>, TV) is positive in all the bands Bn,u for all values of TV and n. Hence, Eb given in (22) has the same 
sign as xsin(TV</>). Following arguments similar to that of the momentum, we find that the binding energy is positive 
in the left part of each band, negative in the right part, and zero at the point <j> = 4>N,n- 

We thus see that for (f> > 0, the momentum and the binding energy are both positive in the left part of each 
band, and they are both negative in the right part. [If <p < 0, we can similarly show that bound states with positive 
(negative) values of P/<j> have positive (negative) binding energy]. In Fig. 2, we show the binding energy Eb as a 
function of (J)/tt for three different values of TV. (We have set fi x 2 = 1 m the figure). We see that Eb is indeed 
positive (negative) in the left (right) part of each band. 

We will now use the technique of continued fractions to study the inverse problem of determining the values of 
TV for which TV-body bound states exist for a given value of <f>. Any positive real number x has a simple continued 
fraction expansion of the form [14] 



no + 



1 



77,1 



(23) 



n 2 + 



where the n^'s are integers satisfying no > 0, and m > 1 for i > 1. The expansion ends at a finite stage with a last 
integer nu if x is rational, and does not end if x is irrational. Given a number x, the integers rii can be found as 
follows. We define xq = x. Then no = [xq], where [y] denotes the integer part of a non-negative number y. We then 
recursively define Xi+i = l/(xi — rij), and obtain rii+i = [xi+i] for i = 0, 1,2, • • ■. If we stop at the fc th stage, we obtain 
a rational number r^ =< no, n\, 112, ■ • • , fife > which is an approximation to the number x. If we write r^ = Pk/qk, 
where p^ and q^ are relatively prime, then it is known that 



Pk_ 
(Ik 



< 



1 



(24) 



for all values of k > 1 [14]. 

Now suppose that we know the expansion 



- = <0,7ii,n 2 ,--- > • (25) 



n 



If we stop at the k stage in this expansion, we obtain pu/qu =< 0, n\, ri2, ■ ■ ■ ,rik >. Eq. (24) then implies that 

\±-*\<\. (26) 

7T qk q% 

We now recall the comment that both the right and the left part of the band B qk iPk have widths which are larger than 
1/qf . Hence Eq. (26) implies that (p/n must lie within the band B qk ^ Pk . We have thus found a value of TV = q^ for 
which an TV-body bound state exists for the given value of <j>. We can generate several such values of TV by stopping 
at different stages k in the expansion given in (25). If 4>/tt is rational, the continued fraction expansion stops at a 
finite stage, so we only obtain a finite number of values of TV in this way. This can also be seen directly from Eq. 
(12). If (p/n = p/q is rational, then at least one of the inequalities in (12) will be violated if TV > q. We thus conclude 
that if <fi/ir is rational, there is only a finite number of values of TV for which a TV-body bound state exists. If (p/ir is 
irrational, then the expansion in (25) does not end, and we can use the procedure described above to find an infinite 
number of possible values of TV for which a TV-body bound state exists. 

To conclude, we have used the ideas of Farey sequences and continued fractions to determine all the allowed ranges 
(bands) of rj for which quantum TV-body bound states exist in the DNLS model. For TV > 3, we find that the TV-body 
bound states can have both positive and negative momentum. Bound states with positive (negative) values of P/rj 
have positive (negative) binding energy. Our work brings the analysis of the quantum bound states in the DNLS 
model to the same level of completion as that of the usual nonlinear Schrodingcr model (where bound states are 
known to exist for all negative values of the coupling constant and all values of TV > 2). 

Bound states with negative binding energy are unusual in the field of intcgrablc quantum models. However, such 
states are known to exist in other areas of quantum physics, such as antibonding states in molecules (see [15] for 
instance). The negative binding energy states that we have found in the DNLS are stable because the model is 
integrable. Presumably, these states would decay if one were to add terms to the Hamiltonian which destroy the 
integrability; any real system would probably have such terms anyway, so it is not clear at the moment if such states 
can be observed experimentally. 

The quantum bound states which exist in the lowest band and have positive binding energy can be related in several 
ways to the solitons which appear in the classical version of the DNLS model which is integrable. A general method 
for relating quantum bound states and classical solitons is described in Ref. [4] . For the case of DNLS model, the 
classical solitons are localized solutions of the equation 



dip f.2 d 2 ip dtp 

ot ax A ox 



* -± = h «I* + ^ M -± , (27) 



with J_ dxvb*vb = TiN . These solitons are known to exist only if < \q\ < w/N [6,13], and they can be obtained by 
taking the 7i > , TV — > 00 limit of the quantum bound states [6] . Another way of relating the classical solitons and 
quantum bound states of DNLS model is indicated in the first paper in Ref. [11]. There it is argued that the classical 
soliton mass, with one-loop quantum corrections, is given by 

M c i - \ [ TV + ?L (TV 3 - TV) ] . (28) 

Comparing this to the mass of the TV-body bound state in Eq. (10), we see that the two agree up to order rj 1 for 
small rj. An interesting problem for future study may be to see if there is a classical version of the bound states in 
the higher bands which we have found in this paper. Since the ranges of values of rj for which these bound states 
exist depend sensitively on TV, going to the limit TV — > 00 with a fixed value of rj may turn out to be a rather subtle 
problem. 
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